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Abstract 

Predictions for light charged hadron production data in the current fragmentation region of 
deeply inelastic scattering from the HI and ZEUS experiments are calculated using perturbative 
Quantum Chromodynamics at next-to-leading order, and using fragmentation functions obtained 
by fitting to similar data from e + e~ reactions. General good agreement is found when the magni- 
tude Q 2 of the hard photon's virtuality is sufficiently large. The discrepancy at low Q and small 
scaled momentum x p is reduced by incorporating mass effects of the detected hadron. By perform- 
ing quark tagging, the contributions to the overall fragmentation from the various quark flavours 
in the ep reactions are studied and compared to the contributions in e + e~ reactions. The yields of 
the various hadron species are also calculated. 
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I. INTRODUCTION 



Due to their high accuracy, data for single hadron inclusive production in high energy 
e + e~ reactions have been used within the framework of the factorization theorem of Quantum 
Chromodynamics (QCD) at leading twist and at next-to-leading order (NLO) to constrain 
fragmentation functions (FFs) for charge summed light charged hadrons (7r ± , and p/p) in 
Refs. [l[ 3, Ij]. The benefits of such an extraction are twofold. First, a test of perturbative 
QCD is provided and consequently imposes a constraint on the strong coupling constant 
a s (Mz) at the Z boson mass scale Mz- Second, since the universality principle of the 
factorization theorem implies that the FFs are independent of the initial state, FFs extracted 
in this way can be used to make predictions for other hadron production processes such as 
those arising from ep reactions in the current fragmentation region and from pp and pp 
reactions. 

Tests of universality were performed in Ref. jlj] by confronting predictions obtained from 
the KKP FF set |l| with corresponding measurements of rapidity (y) and transverse mo- 
mentum (pt) distributions for unidentified light charged hadron production in pp reactions 
at UAl, UA2 and CDF, 7^ reactions at HI and ZEUS, and 77 reactions at OPAL. Within 
the theoretical and experimental errors the description of all data sets was good. However, 
the predictions for the pp reactions have large theoretical errors arising from scale varia- 
tions, and the experimental errors are largest at large pt where additional non-perturbative 
information such as higher twist is expected to be least important. The 7p and 77 reaction 
data suffer from similar problems, but in addition the predictions gain large errors from the 
rather badly constrained parton distribution functions (PDFs) of the photon. Distributions 
in y generally have even larger theoretical errors. More stringent tests of universality were 
performed in Refs. @, Q, [s] through analysis of pseudorapidity and pt distributions from HI 
for the process ep — > e + tt° + X, and in Ref. [6j through analysis of pt distributions from 
ZEUS for the process ep — > e + h + X, which did not require the use of photon PDFs (except 
in the low Q region |9||). The disagreement found with the ZEUS data was reduced in Ref. 
jlo| through resummation of multiple parton radiation at low px- 

In this paper we confront predictions of normalized light charged hadron scaled mo- 
mentum (x p ) distributions with single hadron inclusive production measurements in deeply 
inelastic scattering at the HI and ZEUS 12j, ll3| experiments (the more recent ZEUS 
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data of Ref. 



14| are unfortunately unavailable) at high Q in the current fragmentation re- 



gion, where the detected hadron originates from the fragmentation of a parton at high scale. 
These hadrons can be reliably distinguished from those in the target fragmentation region by 
working in the Breit frame, where the struck quark, which subsequently fragments, moves 
in the opposite direction to the proton remnants, so that x p distributions in the current 
fragmentation region are closely related to x p distributions in any one of the two event 
hemispheres of e + e~ reactions. Consequently, comparison of predictions for ep reaction data 
using FFs constrained from e + e~ reaction data allows for more direct tests of universality. 
Since the data are normalized, uncertainties from the proton PDFs and their perturbative 
evolution are reduced, as well as the dependence on Bjorken x. 

The charge- squared weighted FFs are weighted equally in e + e~ reactions. In particular, 
this implies that FFs for massless d and s quarks cannot be separated if they are not 
separately tagged, so that, since no individually tagged light quark flavour data was used 
in the analyses of Refs. jll, Q, S], additional theoretical constraints on the d quark had to 
be imposed. However, calculation of hadron production processes from proton initiated 
reactions at facilities such as HERA (e + p), the Tevatron (pp), RHIC and the LHC (pp), 
where the charge-squared weighted FFs for quarks of each flavour have an independent 
weighting provided by the PDFs, may demand some degree of knowledge of the individual 
quark FFs, particularly in the light quark sector. 

In the determination of the AKK FF set for light charged hadrons [4], a more phenomeno- 
logical separation of the light quark flavour FFs was pursued using the individually quark 
flavour tagged probabilities measured by the OPAL collaboration 15[. These probabilities 
were constrained by single and double hadron inclusive production measurements for which 
light quarks are favoured, together with the well justified theoretical assumptions of SU(2) 
isospin invariance between u and d quarks for the quark compositions of tt ± , and the branch- 
ing ratios of the Z boson into quark-antiquark pairs of each quark flavour from perturbative 
QCD. Small x p subtleties in double hadron inclusive production are relatively unimportant 
since the data are in the range x p > 0.2. Such a separation should make little difference 
to the current knowledge of FFs, since the SU(2) isospin relation was also used in the 
extraction of the KKP FF set to constrain the d quark. Thus, predictions for ^ data should 
not depend too much on the choice of FF set. The same applies to unidentified light charged 
hadron data, albeit to a slightly less degree, since ^ dominates the sample on account of its 
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low mass. However, the anticipated strange quark suppression in production observed 
in the OPAL experiment gave more realistic production FFs for d and s quarks in the 
AKK set than those in the KKP set, where the FF for s was set equal to the FF for u for 
simplicity. Finally, the light quark separation of the AKK FFs for p/p production may also 
be significant, although it was limited by the large experimental uncertainties of the OPAL 
tagging probabilities. 

These expectations are found to some degree in the comparisons of theoretical predictions 



with pp initiated single hadron inclusive production data [lfj, |l7j from the STAR collabora- 
tion. In Ref. [l6], both the AKK and KKP FF sets lead to similar and good descriptions of 
the it^ yield, while the theoretical prediction gives better agreement (at scale p = p?) with 
the measured p/p production when the p/p FFs are employed from the AKK set than from 
the KKP set. The AKK set for K ± and K s [3] also resulted in an improvement [3] in the 
theoretical description of the K s production measurements of Ref. 17]. 

The paper is organized as follows. We first present the formalism behind our calculations 
in section [III We define the observable we are studying, and give the form of the cross section 
in terms of the FFs to underline the similarities among, and differences between, single 
hadron inclusive production in e + e~ and ep reactions. Then we discuss the modification 
to the cross section when the detected hadron's mass is not negligible, since this effect is 
important at sufficiently small x p and low Q. Section [TTTI contains our comparisons with the 
data, and we examine the uncertainties arising from the arbitrary choice of scale, of PDF set 
and of FF set, as well as the importance of gluon fragmentation and of the detected hadron 
mass effect. Furthermore, although the corresponding data is absent, the contributions 
from the individual fragmenting parton and detected hadron species to the cross section 
are calculated to further determine differences and similarities of the FF sets. In section 
ITVl we present our conclusions. Finally, the appendix gives details on the cuts used in the 
experiments. 



II. THEORETICAL FORMALISM 



We are concerned with the process ep — > e + h + X, where h is a detected hadron and 
X is the remaining unobserved part of the final state, whose kinematic variables will be 
assigned according to the external particles of the general graph in Fig. [TJ The kinematic 
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FIG. 1: General graph for the leading twist contributions to the process e(k)p(P) — > e(k') + h(j>f l ) + 
X. Parallel trios of lines signify unobserved final states. 

degrees of freedom are chosen to be the centre-of-mass (cm.) energy yfs of the initial state 
electron-proton system, which is given by s = (P + k) 2 and which is kept fixed in the 
experiments, the magnitude of the hard photon's virtuality Q 2 = —q 2 , the Bjorken scaling 
variable x = Q 2 /(2P ■ q) and the scaled detected hadron momentum x p = 2ph ■ q/q 2 - The 
normalized cross section (with the s dependence omitted for brevity) takes the form 

f f dQ 2 dx f XpB dx p dOPi ; ton h (x,x p ,Q 2 ) 

^proton h ( . \ _ Jcuts ^ Jx pA V dx v \ V ^ ' ( 

t (^cuts, XpA, %pB) J cnts dQ^dxO^(x,Q 2 ) ' [ ) 

where, for convenience later, we use the shorthand O for d 2 a/(dxdQ 2 ), where "cuts" refers to 
a specified region in the (x, Q 2 ) plane (see the appendix for the various cuts used by HI and 
ZEUS), and where x p a(b) is the lower (upper) edge of the x p bin. The cross section and the 
kinematic variables are frame invariant, and are measured in the Breit frame, defined to be 
the frame where the photon energy vanishes. In this frame the target fragmentation region 
(x p < 0) contains the proton remnants, while the struck parton fragments into the current 
fragmentation region (x p > 0), and the latter process is equivalent to the fragmentation of 
a parton into an event hemisphere in e + e~ reactions. 



A. QCD factorization theorem 

The factorization theorem dictates that the leading twist component of the factorized 
cross section in the numerator of Eq. (0Q) is calculated from processes of the form shown in 
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Fig. [T] and takes the form 

rf £)proton h rl dy rl dz d Qij / Q2 



xx Q 2 ) = f *» f d -lT—(vz^a^)) 

p IJ 



X proton I X 2 \ X P T^h f X P J 
x —Ji "i/ 1 I -^j I iP 



(2) 



In this framework, the incoming parton z has momentum p = (x/y)P and the outgoing 
parton j has momentum p' = (z/x p )ph- ff Toon is the PDF of parton i in the proton, Dj 
is the FF of parton j to the hadron h, % i is the equivalent factorized partonic observable 



given to NLO in Ref. [20|, \i is the factorization / renormalization scale which distinguishes 
the soft from the hard subprocesses and a s (fj, 2 ) = a s (p,)/(2ir). The more commonly written 
form of Eq. (j2J) can be obtained by changing the integration variables to x = x/y and 
x p = Xp/z, which accentuates the role of the PDFs and FFs as probability densities. Using 
the momentum sum rule 

V / dzzD*(z,iJ 2 ) = l, (3) 
h Jo 

the integration over x p from to 1 and the sum over h of Eq. (j2J) yields the factorized cross 
section in the denominator of Eq. ([TJ), viz. 

o*°*»{x, q 2 ) = f ^ & U % -fr ton (-, v 2 ) ■ (4) 

Jx y i v ^ / v \y j 

B. Comparison with e + e~ — » h + X 

We now perform a pedagogical study of the FF dependence of _FP roton h m Therefore, and 
in this subsection only, we work to leading order (LO) (however, all calculations used for 
our numerical analysis of section III will be performed to NLO), where 

d ° 13 f y,z,%M^)) =^(Q 2 )EM^ 7 (Q 2 )5(l-y)5(l-^) (5) 



dz \ /i 2 



and 



& (y, ^, a s {^ = |^(Q 2 ) £ S«el(Q 2 )S(l - y). (6) 

In these expressions, o"o is the cross section for the elastic process e/i — ► e/x for one photon 
exchange in the i-channel, and / indexes the quark of effective electroweak charge e qi (Q 2 ). 
As a further simplification, we will neglect the bin width in Q since F proton h is approximately 
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independent of Q up to 0(1/ In Q) corrections. (However, we will not neglect the bin widths 
in Q in section HTT1 ) The result is 

(cut S ,x, 4 , %B ) E.eiW'^W 2 ) ' (?) 

where Gi{Q 2 ) = J cuts dx xff roton (x,Q 2 ). If the G/ are independent of J, the numerator of 
Eq. ([7]) is equal to the equivalent LO result for e + e~ —> h + X. It is therefore essentially the 
differences between the Gi which distinguishes the two types of observables. The relative 
sizes of the components of the ep cross section where the quark directly connected to the 
electroweak boson is tagged help to determine the relative importance of the fragmentation 
of individual quark flavours in the untagged cross section in ep reactions, and consequently 
to what extent these data could complement the untagged and tagged data from e + e _ 
reactions in understanding fragmentation from the various quark flavours. The ep cross 
section for which quark / is tagged can be obtained by setting the remaining quark charges 
to zero in the calculation, implying that it is scale independent and is given at LO by 
^qjiQ^GiiQ^XpD^Xp.Q 2 ). Then the largest component is the tagged cross section for 
which I = u, due to the valence structure of the proton, the larger charge of the u quark 
relative to the d quark and, to some extent, because the u quark is the most favoured one 
in the production of light charged hadrons. By the nature of data from e + e~ reactions, 
the u quark fragmentation is currently also the most constrained component, in particular 
for tt^ production which constitutes most of the sample, while the most unconstrained 
component is the difference between the d and s quark fragmentations due to their similar 
effective electroweak charges. Consequently, in the absence of sufficiently precise data for ep 
reactions, uds (or equivalently c and b) tagging would therefore be valuable since together 
with e + e~ reaction data it would provide some constraint on the separation between d and 
s quark fragmentation, by virtue of the difference between Gd and G s . So far only charm 



quark tagging through D*^ production measurements 2l| has been performed in the fully 
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inclusive case to obtain ■ Separate u, d and s tagging as performed in e + e~ reactions 
would constitute a further improvement, but may not be possible at present. 

Gluon fragmentation is not so well constrained by e + e~ reactions since the gluon does 
not couple directly to the electroweak boson. Although the proton is an abundant source 
of gluons, this uncertainty is unlikely to contaminate the measurements from ep reactions 
for the same reason. This contrasts with pp and pp reactions, where gluon fragmentation is 



very important because a gluon or quark from one (anti)proton can probe a gluon from the 
other (anti)proton directly and with a much stronger coupling. 



C. Detected hadron mass effect 

The production rate of the detected hadron falls as its mass increases due to the 
reduction in the size of the available phase space. This effect is particularly pronounced 
at small x p and low Q, where mh cannot be neglected relative to the hadron's spatial 
momentum. Treatment of the hadron mass effect in the timelike case was covered in Ref. 



221 ] ; here we derive the modification to Eq. in the spacelike case. The result is essentially 



22J after making the replacement s — > Q 2 , where ^Js is the cm. 



equivalent to that of Ref. 
energy of the e + e~ system. 

In general, the scaling variables of the factorization theorem are given by ratios of the 
light cone momenta. To find the general relation between the true scaling variable of frag- 
mentation and the measured variable x p in the presence of hadron mass, we work in the 
class of frames in which the spatial momenta of the virtual photon and the detected hadron 
are parallel, but is otherwise completely general. It contains, but is not limited to, the Breit 
frame, which is achieved by a boost in the direction of the two momenta. The 3-axis is 
chosen to be aligned anti-parallel with this direction, with no loss of generality. In light cone 
coordinates V = (V+, V~, V T ), where V ± = (l/y/2)(V° ± l^ 3 ) and V T = (V u V 2 ), we then 
have 

Q 2 



q= {-w- q ^ )- (8) 

and the momentum of the detected hadron with non-zero mass in terms of the scaling 
variable £ p = p^/q^, which is invariant with respect to boosts along the 3-axis, is 

This immediately implies that £ p is related to the measured variable x p through 
In Eq. (j2J), the partonic momentum must be chosen as 



p' = -(0,p",0), (11) 



p 



and £ p must replace x p everywhere. The left hand side, dO proton /d£ p , is related to the 
experimentally measured quantity dO pmton h /dx p by 

dO pTOton h , 1 dO h , _ , , ^ 

[x,x p ,Q ) = — 2 — ——(x,£ p {x p ),Q ). (12) 

ax p ]_ _|_ m h ac, p 

Q £,p( x p) 

This normalization of the theoretical cross section agrees with one which has already 

n n 

been proposed [23], and applied in analyses of experimental data [U, |24|], up to terms 
ofO(K/(^ 2 Q 2 )) 2 ). 

In principle, the effect of the initial state proton mass, which is most important at large 
x and low Q, should also be accounted for. However, since the data we will study are mostly 
extracted at small x values, and since this effect modifies the numerator and denominator 
of Eq. (TjQ) in similar ways, we will neglect it. 

III. COMPARISONS WITH HERA DATA 

In this section we present our numerical results for the single hadron inclusive production 
measurements from HI and ZEUS. The kinematic regions of these data are discussed in the 
appendix. In FF fits, uncertainties at small x p , such as higher twist effects, quark and hadron 
mass effects and unresummed soft gluon logarithms in the evolution of the FFs, render the 
theoretical calculations for hadron production data from e + e~ reactions unreliable when the 
scaled momentum, given in this case by x p = Ipnl y/s, where ph is the cm. momentum of the 
detected hadron, falls below 0.1. Because of the resulting uncertainties in the FFs at small 
x p , and because ep reaction data suffer from similar uncertainties at small x p , we only study 
ep reaction data for which x p > 0.1. Cross sections are calculated to NLO in the MS scheme 
using the CYCLOPS software [25]. We set the number of active quark flavours rif — 5. 



To account for the initial state proton, we use the CTEQ6M PDF set of Ref. 26j unless 

(5) 

otherwise stated. We use their value Aq^, d = 226 MeV. Although this does not coincide 
with the values at which the various FF sets are obtained, within this range of values the 
dependence on Aq^, d is rather small. The factorization / renormalization scale is chosen 
as fi = Q unless stated otherwise. The detected hadron's mass is set to zero unless 
otherwise stated. 
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A. Scaled momentum distributions 



In this subsection we compare theoretical predictions with single hadron inclusive pro- 
duction x p distributions measured by HI [11] (see Fig. [221 for the kinematical constraints) 



and ZEUS 



121 ] (see Fig. |2"B"|) . The predictions generally agree well with the ZEUS data (Fig. 



[2]). The predictions using the Kretzer FF set 2[ are similar to those in Ref. )lo| . where 
the CTEQ5M1 PDF set was used. A similar comparison was performed in Ref. [27| using 



the BKK FF set 



and the agreements were good when the CTEQ3M and MRSA' PDF 
sets were used. For both the HI (Figs. [3] and Hj) and ZEUS data, the predictions using 
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FIG. 2: Comparisons of theoretical predictions using the AKK, Kretzer and KKP FF sets with 
the x p distribution from ZEUS \\" 



the KKP FF set are the most gradual in x p , while the Kretzer predictions are the steepest. 
The predictions from the AKK and Kretzer sets are quite similar, particularly at large x p 
and for all x p values of the high Q HI data (Fig. H]). The uncertainty from the freedom 
in the choice of FF set is largest at large x p , since the data from e + e~ reactions is most 
inaccurate and most scarce at large x p . The predictions for the low Q HI data (Fig. [3]) show 
an undershoot at large x p . This behaviour may result from unresummed logarithms at large 
x p in the partonic cross section, since resummation tends to enhance the cross section. The 
overshoot from the low Q HI data at small x p may be due to the theoretical errors in ep 
reaction data discussed above. Indeed, better agreement is found at small x p with the high 
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Q HI data (Fig. where resummation is less necessary and where higher twist and mass 
effects are significantly reduced. 

We now study various modifications to the predictions for the low and high Q HI data in 
order to understand the effect of increasing Q on the theoretical and propagated experimental 
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FIG. 3: As in Fig. El for the low Q HI x p distribution [11] 
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FIG. 4: As in Fig. O for the high Q HI x p distribution 
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errors. First we modify our theoretical approach to incorporate the detected hadron mass 
according to the method of subsection III CI Since the hadron sample is dominated by pions, 
the "average" hadron mass is expected to be around = 0.2 — 0.3 GeV. However, to 
exaggerate the effect of hadron mass for illustration, we choose the larger value rrih = 0.5 
GeV. At small x p , this effect improves the description of the low Q HI data (Fig. UJ), while 
making negligible difference to the high Q HI data (Fig. [6]) over the whole x p range. However, 
this improvement should not be taken too seriously, since other low Q, small x v effects may 
also be relevant. In addition, the FFs from the various sets are artificially suppressed at 
small x p since the hadron mass effect was not accounted for in the analyses of Refs. jl, 2, J]. 




FIG. 5: As in Fig. El using only the AKK FF set. The modifications to the default predictions 
(solid line) arising from the replacement of the CTEQ6M PDF set by the MRST2001 PDF set of 
Ref. [29|, from the removal of the evolved gluon, and from the incorporation of the hadron mass 
effect are shown. 



The error due to the freedom in the choice of PDF set, which we determine by calcu- 
lating the predictions using the MRST2001 PDF set 29|, is rather small, particularly for 
intermediate x p values and for the high Q data. 

The gluon contribution (also shown in Fig. [5]) is clearly negative, although the evolved 
gluon FF is positive. This quantity is calculated by setting the evolved quark FFs to zero 
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FIG. 6: As in Fig. [5J for the high Q HI x p distribution. 

in ^P roton h . Although it is scheme and scale dependent, its definition is the same for all 
3 FF sets, and therefore its variation with respect to the choice of set is due only to the 
different choices of e + e~ reaction data used in the fits. In general, the gluon fragmentation 
is unimportant, particularly away from the smaller x p range and for the high Q HI mea- 
surements. For the low Q HI data (Fig. [7]), the uncertainty from the gluon fragmentation 
from its average is about ±4% at the smallest x p range and about ±2% at x p ~ 0.5. This 
reduces to ±2% and ±1% at the same respective x p values for the high Q HI data. The 
gluon FF is least important for the Kretzer predictions, and most important for the AKK 
ones. 

For the low Q HI data, the uncertainty from the freedom in the scale choice is largest 
at the smaller and larger x p values (Fig. [8]). In addition, since Q is low, our neglect of 
charm quark threshold effects is expected to contribute significant errors at small x p . These 
uncertainties are most likely dominated by unresummed logarithms at small and large x p 
discussed above since, for the high Q HI data (Fig. [9]), the error at small x p is much smaller, 
while the error at large x p still remains sizeable. In general, increasing the scale steepens 
the drop in the cross section with increasing x p . 

To determine how the relative importances of the fragmentations from the various quark 
flavours differ between the ep and e + e~ reaction data, we study the quark flavour tagged 
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FIG. 7: Ratios of the evolved gluon contribution to the HI low Q measurement to the full mea- 
surement, calculated using the AKK, Kretzer and KKP FF sets. 
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FIG. 8: As in Fig. [5l for the modifications arising from scale variation. 



components of the cross section in Fig. [TU] (the low Q predictions are not considered since, 
as we have just seen, the theoretical errors are larger), and compare with the quark tagged 
cross sections versus x p in e + e~ reactions (Fig. [IT]) . The latter reaction was calculated 



14 



using the method of Ref. J4]. As anticipated from the valence structure of the initial state 
proton in subsection III Bl the contribution to the overall fragmentation from the u quark 
fragmentation constitutes a significant amount (50% or more) of the HI and ZEUS data, 
while the contribution from d quark fragmentation is much less. In the e + e~ reaction data, 
the u and d quark fragmentations feature in roughly equal proportions since their FFs and 
electroweak charges are similar, and together contribute 50% or less to the production. 
Fragmentation from s quarks is more important in e + e~ reactions, particularly at large 
x p , while fragmentation from the c quark constitutes similar fractions in both reactions. 
Generally, the part of the fragmentation arising from the b quark is small due to its small 
charge and high mass. However, while it can be relevant in e + e~ reactions at small x p , 
it is always negligible in ep reactions due to its low density in the proton. In the single 
hadron inclusive production data for ep reactions at large x p , the contribution to the overall 
fragmentation from the s quark is more important in the predictions of AKK and Kretzer 
than in the KKP predictions. On the other hand, as expected, all FF sets lead to similar 
contributions from the u, d + s and c quark fragmentations, and b quark fragmentation is 
always negligible. 

The relative importances of the fragmentations into the various light charged hadrons in 
ep reactions can be determined from the composition of the detected hadron sample with 
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FIG. 9: As in Fig. [6l for the modifications arising from scale variation. 
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respect to the hadron species (Fig. [T2|) . (The Kretzer p/p FFs were calculated by subtraction 
of the tt^ and K FFs from the FFs for all light charged hadrons — no p/p production data 
was used in the extraction of the Kretzer FFs.) The uncertainty in the different yields is 
estimated by the spread of the results for the different FF sets, and is largest at large x p 
and smallest at intermediate x p . The AKK and Kretzer sets give rather similar descriptions 
of the ir^ and yields for all x p values shown, while the KKP set gives larger yields at 
large x p . Fragmentation to p/p at large x p , where all three predictions differ considerably, 
is clearly difficult to calculate reliably. The yield in ep reactions (Fig. [T3j) . as for 
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FIG. 10: The ratios of the quark tagged components of the cross section to the untagged cross 
section for the high Q HI data, using the AKK, Kretzer and KKP FF sets. The lowest 3 curves 
show the contribution from the u quark tagged component only, the next 3 curves above the sum 
of the u and d components, the next 3 u, d and s etc. 

e + e~ reactions (Fig. fT4l) . is the largest one due to the low mass of the charged pion. The 
fraction of is slightly larger in e + e~ reactions than in ep reactions, possibly because s 
quark fragmentation is more important in the former data: The most important source of 
K ± is the s quark, since the other favoured quark, u, has to extract a heavier s quark from 
the sea. 
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FIG. 11: The ratios of the quark tagged components of the e + e — ► h + X cross section, where h 
is any light charged hadron, to the untagged cross section, at y/s = 91.2 GeV and using the AKK 
FF set. 
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FIG. 12: The individual hadron species constituting the sample for the high Q HI data, using the 
AKK, Kretzer and KKP FF sets. Each curve is for a single hadron, not a summation of hadrons. 
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FIG. 13: The ratios of the individual hadron species constituting the sample for the high Q HI 
data, using the AKK, Kretzer and KKP FF sets. 




FIG. 14: The ratios of the individual hadron species constituting the sample for the x p distribution 
of the e + e~ — > h + X cross section, where h is any light charged hadron, to the cross section for 
the full sample, at yfs = 91.2 GeV and using the AKK FF set. 
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B. Distributions in photon virtuality 



Next we compare theoretical predictions with the single hadron inc 



measurements at various Q values from HI 



ll| (see Fig. |22} and ZEUS 



usive production 
see Fig. 



The predictions agree well with the ZEUS data (Fig. IT51) . except for, at low Q, the overshoot 
at small x p and the undershoot at large x p . Similar behaviour is found with the less precise 
HI data (Fig. [T6|) . Note that the theoretical predictions are rather constant over the whole 
Q range of both data sets, as foreseen in subsection III B[ Except at the lower Q and smaller 
x p region, the AKK predictions tend to be closer to the Kretzer predictions than to the 
KKP ones. 

The hadron mass effect brings the prediction closer to the data (Figs. [T7] and [T51) at low 
Q and small x p . In this region charm threshold effects are expected to be important, and 
this may explain the large average hadron mass required to obtain convergence of the theory 



24| by essentially 



with the data. Good agreement with the HI data was obtained in Ref. 
choosing nih = 0.66 GeV. 

The uncertainty from the freedom in the choice of PDF set for the proton (Fig. [T5j) is 
everywhere insignificant. At smaller x p values, the gluon fragmentation and the uncertainty 
with respect to the arbitrary scale choice (Fig. [T9I) become less relevant with increasing Q, 
and are unimportant for all Q at the other x p values. The large deviation of the prediction 
for fi = Q/2 (dotted line) from the one for // = Q is caused by the vanishing of the c quark 
FF below threshold. This behaviour is not physical since we have neglected charm mass 



effects. The procedure for incorporating these effects is given in Ref. [30], which amounts 
to retaining the heavy quark mass dependence in the heavy quark flavour creation from 
photon-gluon fusion, and using the same scaling variable that results in the latter process 
for the heavy quark flavour excitation. Furthermore, the matching conditions of Ref. 31] 
must be imposed on the FFs at the quark flavour thresholds. In any case, our results at low 
Q suffer other theoretical errors mentioned earlier, such as higher twist. 

At the lower Q values, the second most important source of fragmentation after the u 
quark is the fragmentation from the c quark (Fig. I2UI) . although, for the HI data, this falls 
with rising Q until the d quark fragmentation becomes more important. 

The relative yield of each hadron species does not change significantly with Q (Fig. I2T1) . 
which is expected from perturbation theory at high Q. 
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FIG. 15: Comparisons of theoretical predictions using the AKK, Kretzer and KKP FF sets with 
the ZEUS data [2]. Each data set is measured in a specific x-bin and, together with its predictions, 



is shifted upwards relative to the one below by the indicated value for A. 



IV. CONCLUSIONS 



We have performed a comprehensive analysis of single hadron inclusive production data at 
HERA, by calculating the theoretical predictions using FF sets that were obtained by fitting 
to accurate e + e~ data. In general, good agreement was found using the AKK, Kretzer and 
KKP FF sets. However, at low Q and small x p the predictions overshoot the data, a problem 
which is partially remedied by including the detected hadron mass effect. Unresummed soft 
gluon logarithms may also contribute to this discrepancy, as suggested by the increasing 
variation with respect to the scale for decreasing x p and Q, as well as higher twist and quark 
mass effects. A more complete treatment which takes into account all these effects is needed 
to improve the understanding of fragmentation in this region. At large x p , an undershoot 
occurs in the HI data, which may be avoidable by resumming the logarithms at large x p , 
as suggested by the small rise in the variation with respect to the scale as x p increases and 
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FIG. 16: As in Fig. [J5l for the HI data 



Q decreases. As is the case for e + e~ reactions, gluon fragmentation is not important in ep 
reaction data, at least for sufficiently large x p and Q. 

Fragmentation from the u quark gives the largest contribution to the overall fragmentation 
in ep reactions, followed by c and d quark fragmentation (c quark fragmentation being more 
important at lower Q), and finally from the s quark. Fragmentation from the b quark is 
negligible. This should be contrasted with the situation in e + e~ reactions, where the u and d 
quark fragmentations are of similar importance, while fragmentation from c is less important 
than from s, and b quark fragmentation has some relevance at smaller x p . Therefore, even 
sufficiently accurate untagged e + e~ and ep reaction data taken over a large enough range of 
the kinematic variables would improve the constraints on the individual quark flavour FFs, 
although quark tagging is more valuable for this purpose. 

The fractional yields of each of the light charged hadron species in ep reactions depend 
strongly on x p , but to a much lesser extent on Q. They are similar to the fractional yields 
in e + e~ reactions. 

Relative to the experimental accuracy of the data sets, the AKK and Kretzer predictions, 
as well as their quark tagged components and n ± and K ± yields (but not the evolved gluon 
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FIG. 17: As in Fig. [16] for the data measured in the interval 0.1 < x p < 0.2, using only the AKK 
FF set and for different values of rah- 

fractions and the p/p yields) are very similar for all data considered. At the time of writing, 
the HI and ZEUS collaborations are planning an extraction of very accurate data using, 
respectively, improved triggering and higher luminosity, which could allow for a comparison 
of the reliablility of the FF sets. 

APPENDIX: EXPERIMENTAL CUTS 

In this section we present the regions in (x, Q 2 ) used by the HI and ZEUS collaborations 
from which the measured cross sections are extracted. These regions are bounded according 
to cuts on x, Q 2 , the squared cm. energy of the virtual photon-proton system, 



and the fraction of the energy of the initial electron (we do not distinguish between electrons 
and positrons) which is lost in the rest frame of the proton, 




(A.l) 



P-q Q 2 



(A.2) 



y = 



P ■ k xs 
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FIG. 18: As in Fig. [T5| using the AKK FF set. The modifications to the default predictions (solid 
ine) arising from the replacement of the CTEQ6M PDF set by the MRST2001 PDF set of Ref. 
29 ], from the removal of the evolved gluon, and from the incorporation of hadron mass effects are 

shown. 



A lower bound on the scattered electron's energy 

where E is the energy of the initial electron, is sometimes imposed to prevent the scattered 
electron being falsely identified with isolated low energy deposits in the calorimeter while 
the true scattered electr on p asses undetected down the beam pipe. The HI collaboration 



imposes additional cuts 



32j on the angle of deflection of the electron and struck parton, 



respectively 9 e and 6 P , to maintain good detector acceptance. In the laboratory frame, these 
are given in terms of x and Q 2 by 

_ xs (AE 2 - Q 2 ) - AE 2 Q 2 
C ° S e ~ xs (4£2 + Q2) - 4 £2Q2 
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FIG. 19: As in Fig. [TBI for the modifications arising from scale variation. 



and 



xs(xs - Q 2 ) - AE 2 Q 2 



cos0„ = — ) ^ n ^L ■ (A. 5) 

p xs(xs - Q 2 ) + AE 2 Q 2 v ; 

The cuts used by HI in Ref. UJ are shown in Fig. [221 The lower bound on 8 e = 10° does 



not bound any of the regions of measurement. The Q bins in this analysis are very narrow 
and are not shown. The cuts used by ZEUS in Ref. [T2J are shown in Fig. [23j It is clear 
that only the cuts in W and Q border the region of measurement. Figure EH shows that the 
upper bound on y is irrelevant in the extraction of the data of Ref. 13]. 
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FIG. 24: Cuts in the (x,Q 2 ) plane used in the ZEUS analysis of Ref. [13]. The measured regions 
are indicated by the labels "ZEUS". E = 27.5 GeV and y/s = 300.3 GeV. 
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